SCIENTIFIC 

REPpRTS 



»H5 




SUBJECT AREAS: 
ATTOSECOND SCIENCE 
NONLINEAR OPTICS 
APPLIED PHYSICS 
X-RAYS 



Efficient high-order harmonic generation 
boosted by below-threshold harmonics 

F. Brizuela 1 , C. M. Heyl\ P. Rudawski 1 , D. Kroon 1 , L Rading 1 , J. M. Dahlstrom 2 , J. Mauritsson 1 , 
P. Johnsson', C. L. Arnold 1 & A. L'Huillier' 

department of Physics, Lund University, P. O. Box 1 1 8, SE-22100 Lund, Sweden, 2 Atomic Physics, Fysikum, Stockholm University, 
AlbaNova University Center, SE-10691 Stockholm, Sweden. 



Received 

12 February 2013 

Accepted 

13 February 2013 

Published 
11 March 2013 



Correspondence and 
requests for materials 
should be addressed to 
A.L. (anne.lhuillier@ 
fysik.lth.se) 



High-order harmonic generation (HHG) in gases has been established as an important technique for the 
generation of coherent extreme ultraviolet (XUV) pulses at ultrashort time scales. Its main drawback, 
however, is the low conversion efficiency, setting limits for many applications, such as ultrafast coherent 
imaging, nonlinear processes in the XUV range, or seeded free electron lasers. Here we introduce a novel 
scheme based on using below-threshold harmonics, generated in a "seeding cell", to boost the HHG process 
in a "generation cell", placed further downstream in the focused laser beam. By modifying the fundamental 
driving field, these low-order harmonics alter the ionization step of the nonlinear HHG process. Our 
dual-cell scheme enhances the conversion efficiency of HHG, opening the path for the realization of robust 
intense attosecond XUV sources. 



The interaction of intense laser pulses with atomic or molecular gas media leads to the generation of 
harmonics of the laser light, up to very high orders 1 . These harmonics are locked in phase, giving rise to 
attosecond bursts of XUV light. The simplicity of the experimental technique, together with the progress in 
ultrafast laser technology, has promoted HHG sources as essential tools in many laboratories; opening, in 
particular, the field of attosecond science 2 . However, HHG suffers from low conversion efficiency, owing partly 
to phase mismatches in the nonlinear medium that prevent efficient build up of the macroscopic field 3-6 , but 
mostly to the weak response of the individual atoms to the field. 

The atomic response to an external driving field can be described by a three-step model [Fig. 1(a)]: First, a 
bound electron tunnel-ionizes into the continuum; second, it is accelerated by the laser field; and finally, it 
recombines with the parent ion upon field reversal, emitting an XUV photon 7,8 . The electron trajectories can 
be grouped in two families, named the long and the short, depending on the excursion time of the electron and 
generated in intervals II and III of Fig. 1(a), respectively. The most interesting from a practical point-of-view are 
the short trajectories, which lead to collimated and spectrally narrow emission. Unfortunately, these trajectories 
start at times close to the zero-crossings of the driving electric field, suffering from very low quantum-tunneling 
probability. 

Altering the driving electric field at the subcycle level 9 provides a way of modifying the single atom response. 
This has been investigated mainly by adding the second harmonic field 10 " 13 , thus breaking the symmetry between 
consecutive half cycles. In contrast, odd-order harmonics modify the HHG process while maintaining the half- 
cycle symmetry. In a pioneering work, Watanabe and coworkers 14 investigated the influence of the third har- 
monic (TH) on single ionization and HHG in Ar, obtaining an enhancement of up to a factor of ten for the 
27-31 harmonics. Also, a few theoretical works discuss the influence of the TH on the enhancement of the 
yield 1516 and/or the extension of the cutoff energy 17 " 19 . Another approach to enhance the signal by modify the 
single atom response is to control the time of ionization by using attosecond pulse trains to initialize the three-step 
process via single photon absorption 20 " 23 . 

In this letter, we demonstrate a simple and robust, yet powerful enhancement scheme based on a dual gas-cell 
setup [Fig. 1(b)] . We study HHG in neon using a high-energy (—20 mj), near-infrared fundamental field, loosely 
focused in a long gas cell, resulting in high-order harmonics in the 40-100 eV range, with a typical energy of 10 nj 
per harmonic order. The addition of a high-pressure Ar gas cell before the generation cell produces a large 
enhancement in the Ne signal, as seen in Fig. 1(c). We experimentally and theoretically show that the observed 
enhancement is due to below-threshold, low- order harmonics which modify the fundamental field in such a way 
that the contribution of the short trajectories is increased. 
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Figure 1 | HHG in a dual gas-cell, (a) Schematic classical trajectories for a sinusoidal driving field (red line). The colors indicate the return energy of the 
electrons in units of the ponderomotive energy U p . Modifying the driving field by adding an odd harmonic field (blue line) can lead to an enhanced 
ionization probability for short trajectories (interval III) while suppressing the ionization of non-contributing electrons (intervals I and II), as indicated 
by the arrows, (b) Schematic experimental setup. Low-order harmonics generated in the seeding cell co-propagate with the fundamental into the 
generation cell and modify the HHG process, (c) Comparison of a typical HHG spectrum from a neon-filled generation cell obtained using only the 
fundamental field; and a spectrum obtained combining the fundamental field with low-order harmonics generated in the argon-filled seeding cell. In the 
latter case, the harmonic yield for the plateau harmonics is enhanced while the cutoff energy and the divergence are reduced. 



Results 

In our experiment, the generation cell is placed approximately at the 
laser focus while the seeding cell is located a few centimeters before 
(see Methods). The gas pressures in the cells can be independently 
adjusted and are typically a few mbar in the generation cell (Ne) and 
up to tens of mbar in the seeding cell (Ar). In Fig. 2(a-c), HHG 
spectra from neon are plotted as a function of the seeding pressure 
for three different driving intensities. When no gas is present in the 
seeding cell, standard Ne spectra are obtained. As the seeding pres- 
sure increases, the signal from the neon cell decreases until it is 
almost completely suppressed. At higher pressures, the neon spectra 
reappear and are significantly enhanced in the 50 - 80 eV region 



while the maximum photon energy slightly shifts to lower harmonic 
orders. 

Figure 2(d) shows harmonics generated in the seeding cell. 
Harmonics with energies above the ionization threshold are not 
present at pressures where the enhancement in the generation cell 
occurs, and therefore are not responsible for the signal boost through 
single-photon ionization 20 " 23 . At these pressures, only low-order har- 
monics are efficiently generated in the seeding cell, indicating that 
they are responsible for the seeding process. 

In order to validate our interpretation, we performed numeri- 
cal simulations for both cells. In the generation cell, we simulated 
the seeded HHG process using the strong-field approximation 




Photon energy (eV) 

Figure 2 | Experimental HHG spectra, (a-c) Spectra from the generation cell as a function of the pressure in the seeding cell at three driving intensities 
2.7, 3.5, and 4.4 X 10 14 W/cm 2 , respectively. The spectra were obtained using argon in the seeding cell and neon, at a fixed pressure, in the generation cell. 
The data were normalized to the most intense enhanced neon spectrum, (d) Low- (3-7) and high-order harmonics from the seeding cell as a function of 
Ar pressure. The dotted lines indicate regions measured independently with different detectors. Each region was normalized to the highest intensity in the 
corresponding spectral range. 
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' (see Methods). The total field can be written as 



E(t)=E 0 



sin(cof)+ r q sin (qcot-\-A(p i 



(i) 



where E 0 is the amplitude of the fundamental field, w its frequency, I p 
the ionization energy, r q the ratio between the fundamental and qth 
harmonic field, and A(p q their relative phase. Although all harmonics 
below the ionization threshold of Ar may influence the enhancement 
phenomenon, we considered only the TH, which is the most intense 
one (we omit the subscript 3 below). A simulated HHG spectrum in 
neon with \r\ 2 = 0.01, is shown in Fig. 3(a) as a function of A(p. A 
relative phase of ~ 1 rad leads to an enhanced ionization probability, 
since the electrical field is increased at the time where the short 
electron trajectories are born [interval III in Fig. 1(a)]. Further- 
more, the electric field amplitude is reduced around the peak of 
the fundamental field leading to suppressed probability for non- 
contributing trajectories (intervals I, II) and to an improved mac- 
roscopic situation since plasma dispersion and depletion effects are 
minimized 4,5 . When A(p ~ 1 ± n, the situation is reversed and HHG 
is suppressed compared to the unseeded case. 



We experimentally confirmed the dependence of the HHG signal 
on A(p by studying HHG using a combination of the fundamental 
and the TH generated in a crystal 14 . To control the delay between the 
two fields, we used a Michelson interferometer with the TH produced 
in one arm. Our results, plotted in Fig. 3(b), show a strong delay 
dependence of the harmonic yield. However, we could not increase 
the overall HHG efficiency compared to the dual-cell scheme, since a 
large fraction of the fundamental field was needed for the TH gen- 
eration and consequently lost for HHG. 

In the seeding cell, we examined the pressure dependence of both 
low-order and high-order harmonic generation. Our calculations 26 
confirm the experimental observation that HHG in Ar peaks at a 
certain pressure (—10 mbar) which corresponds to optimized phase 
matching 27 , while below-threshold harmonics continue to increase 
up to pressures as high as 100 mbar. We also investigated the pro- 
pagation of the fundamental and TH fields in a high pressure cell 28 
(see Methods). This allowed us to examine their phase relation after 
the seeding cell and to eliminate the relatively weak reshaping of the 
fundamental field in our experimental conditions as possible cause 
for the enhancement. As Fig. 3(c) shows, for high enough seeding 
pressures, Acp will be between 0 and 2 radians during part of the laser 
pulse, leading to a gated enhancement mechanism. 
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Figure 3 | Influence of the relative (co, 3co) phase in HHG. (a) SFA spectra 
as a function of A(p in the generation cell, normalized to the unseeded 
spectrum. Only the contribution of the short trajectory is considered. An 
effective grating response is included to mimic the experimental 
conditions, (b) Experimental results with the TH generated in a crystal, 
normalized to the highest signal, (c) Propagation simulations in the 
seeding cell: A(p at the exit of the cell as a function of time for different 
pressures. 



Discussion 

As in any enhancement scheme, a key question is whether our 
method is advantageous over "usual" HHG optimization, which 
can be achieved for example by using looser focusing, optimizing 
the position of the focus in the cell, or adjusting the pressure in the 
gas cell 4,29 ' 30 . Ideally, one would like to compare optimized HHG and 
optimized seeded HHG for a given fundamental pulse energy. This is 
not easy to realize experimentally, so we choose to benchmark seeded 
HHG against optimized unseeded HHG, with —10 nj at 63 eV (41st 
harmonic). 
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Figure 4 | Optimization of HHG. (a) 41st harmonic energy as a function 

of the driving intensity for seeded (red) and unseeded (blue) HHG. 

Unseeded HHG is optimized at the maximum intensity. 

(b, c) Corresponding experimental spectra at 3.5 and 4.4 X 10 14 W/cm 2 , 

respectively. 
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Figure 4(a) compares the 41st harmonic signal in the seeded and 
unseeded cases as a function of the driving intensity. The intensity 
required for saturating seeded HHG is only half that needed for 
unseeded HHG. This explains the reduction of the cutoff energy 
and the lower divergence for the harmonics. The enhancement factor 
depends on the driving intensity [Fig. 4(b,c)]. For the 41st harmonic, 
it varies from five at 3.5 X 10 14 W/cm 2 (and even higher at lower 
intensity) to two at 4.4 X 10 14 W/cm 2 . By further optimizing seeded 
HHG (e.g. by changing the focusing conditions) one should be able 
to obtain an even larger increase compared to unseeded HHG. The 
higher efficiency together with the lower divergence leads to a 
brighter source of XUV light. 

In summary, we have studied the effect of seeding HHG using 
harmonics generated in a separate gas cell and showed that low-order 
harmonics are responsible for the resulting enhancement. The com- 
bined electric field preferentially enhances the short trajectories 
while suppressing depletion and plasma dispersion effects. The 
required phase difference between the fundamental and the low- 
order harmonics is obtained by adjusting the pressure in the seeding 
cell, thus modifying the free-electron dispersion. Our method is not 
limited to the gas combination presented here. Experimentally we 
have observed an increased harmonic yield for a variety of gas com- 
binations, and even when the same gas is used in both cells. Our 
simulations show that the enhancement can be scaled far above one 
order of magnitude by increasing the low- order harmonic intensity, 
for example by using longer cells, higher pressures or gases with 
higher nonlinearities. This also leads to a shorter temporal gate, of 
interest for single attosecond pulse generation. 

Methods 

Experimental setup. The harmonics were generated using 45 fs pulses, centered at 
800 nm. The gas cells used in this setup were 1 cm long with a diameter of 1 mm. The 
injection of gas into the cell was synchronized with the laser repetition rate (10 Hz) 
and the delay between the gas injection and the laser pulse was optimized for each cell. 
In the experiments, seeding cell pressure and pulse energy were the parameters 
investigated. The generation cell pressure was set for the best phase- matching 
conditions for Ne at the highest laser intensity intensity (4.4 X 10 14 W/cm 2 ), 
corresponding to less than 10 mbar. The focus position was adjusted in order to 
optimize HHG in the generation cell. The cell separation was 15 mm with the 
generation- cell located at focus (f = 4 m). Nevertheless, larger separations, up to 
50 mm yielded similar results. The cells were mounted on motor-controlled XYZ 
stages with motorized XY tilt capabilities. The cells could be removed completely 
from the IR field. A CCD camera was used to align each cell to the laser and 
observation of the spectra at the best phase-matching conditions were used to 
evaluate the tilt of each cell. The same Ne spectra from the generation cell could be 
obtained though the evacuated seeding cell or with the seeding cell removed from the 
beam path. The same was true for the seeding cell where Ar spectra could be obtained 
under both conditions. The pressure and intensity controls were automated to scan 
the region of interest. At each experimental condition 10 single- shot spectra were 
measured and averaged. The harmonic orders were calibrated using the absorption 
edge of an Al-foil filter. The fundamental intensity was estimated from the cutoff of 
the unseeded Ne spectra. 

Numerical simulations. Generation cell. The influence of a weak third harmonic field 
on the HHG process was simulated by solving the time dependent Schrodinger 
Equation within the strong field approximation. The quasi- classical action for the 
electron motion in the continuum 




is calculated for a combined vector potential of the fundamental field and a weak 
parallel auxiliary field consistent with the field definition in Eq. (1). t 0 and t 
correspond to the tunneling and recombination times for an electron with canonical 
momentum p, I p is the ionization potential, and A the vector potential of the field. We 
approximate the HHG dipole as 24 

(3) 

xd x (p st (t,T)-A x (t-x))E(t-T) 
x exp[ — iS st (t,T)]F(z) + c.c., 

where a stationary phase approximation is performed over momentum, with 
p st (t, t) = [E(t) — E(t — t)]/t, where x = t — t 0 is the excursion time in the continuum. 



We also insert a filter function F(t) to select the short trajectory: F(t) ~ 1 for 
t < 0.65 T and F(t) ~ 0 for t > 0.65 T, where 0.65 T corresponds to the position of the 
cutoff. The integral in Eq. (3) is then evaluated numerically on a finite grid followed by 
a numerical Fourier transform for the dipole emission. 

Seeding cell. We performed calculations which combine the solution of the 
time-dependent Schrodinger equation in a single-active electron approximation 
and propagation in a partially ionized medium 26,5 using a slowly- varying envelope 
approximation. Our main goal was to examine the influence of the pressure both for 
low-order and high-order harmonic generation in conditions mimicking the 
experiment. We found a maximum for HHG at around 10 mbar, while 
below-threshold, low-order harmonics which are not reabsorbed in the medium 
continue to increase up to very high pressures (100 mbar). 

The generation of the third harmonic in the seeding cell was simulated using a 
(3 + 1) -dimensional, unidirectional, nonlinear envelope equation 28 . The complete 
frequency dependent dispersion relation is considered, enabling to propagate the 
fundamental and the third harmonic simultaneously. It is numerically integrated 
using a split- step technique, where the linear contributions, such as dispersion and 
diffraction are treated in A;-transverse frequency space, while the nonlinear part, 
taking into account the Kerr effect, third-harmonic generation as well as plasma 
dispersion and plasma defocusing is treated in normal space. The method is described 
in detail in 28 . The calculated phase variation is mainly due to plasma dispersion 
effects. There are also small contributions from the geometrical phase acquired along 
the seeding cell as well pressure- dependent third harmonic phase matching. 
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